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Abstract 

We show that the resummation technique developed recently for heavy- 
light systems and the conventional approach based on the factorization of 
Wilson loops are equivalent. The two methods lead to the same results of 
the resummation of large corrections in inclusive B meson decays. 



1 



1. Introduction 



Perturbative QCD (PQCD) including Sudakov effects |], [| has been pro- 
posed to be an appropriate approach to exclusive heavy meson decays in the 
maximal recoil region of final-state hadrons || |], |5|, ||. The progress is 
attributed to the resummation of large radiative corrections in heavy-light 
systems such as a B meson containing a light valence quark, through which 
important higher-order and higher-power contributions are included into fac- 
torization theorems. Recently, the resummation technique was extended to 
the inclusive semileptonic decay B — > X u lv [[/]], whose behavior near the high 
end of charged lepton energy is crucial to the extraction of the Cabibbo- 
Kobayashi-Maskawa matrix elements \V u b\ in the standard model and to the 
detection of new physics. It has been observed that large perturbative cor- 
rections appear at the end point of the spectrum || || [l0| , which must be 
organized to all orders in order to have a reliable PQCD analysis. 

The above approach to the resummation of large corrections in inclusive 
processes differs from the conventional one. In the study of deep inelastic 
scatterings and Drell-Yan productions |TTJ, large soft corrections are factor- 



ized systematically into path-ordered exponentials, or Wilson lines [12, 13, 14 



along the classical trajectories of quarks. Korchemsky and Sterman [L5| ap- 
plied this Wilson-loop formalism to the study of end-point singularities in 
inclusive B meson decays. The Wilson lines are absorbed into a heavy quark 
distribution function S, and give its evolution governed by renormalization 
group (RG) equations. The emitted light quark produces a jet of collinear 
particles at the end point due to its large energy and small invariant mass 

0- 

In this letter we shall show that the resummation technique is equivalent 
to the Wilson-loop formalism, if relevant evolution scales are chosen properly. 
The two methods give the same resummation results for the inclusive decays 
B — > X s 7 and B — > X u lv up to next-to-leading logarithms. Therefore, this 
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work provides a non-trivial check on the two different approaches. 



2. The Radiative Decay B —> X s j 

We apply the resummation technique to the inclusive radiative decay 
B — > X s 7, which occurs through the transition b — > S7 described by an 
effective Hamiltonian |Tj|. The basic factorization of this process is shown 
in fig. la, where the bulb represents the b quark distribution function S(k), 
k being the momentum carried by light partons in the B meson, k 2 m 0. 
k has a plus component k + and small transverse components k^, which 
serve as the infrared cutoff of loop integrals for radiative corrections. The 
momentum of the b quark is then P — k, P = i71b/V2(1, 1,0) being the 
B meson momentum in terms of light-cone components, which satisfies the 
on-shell condition (P — k) 2 « m 2 . Here tub and m& are the B meson mass 
and the b quark mass, respectively. The b quark decays into a real photon 
of momentum q and a s quark of momentum P s , which is regarded as being 
light. Assuming that the nonvanishing component of q is q + , we have P s = 
P-k-q=(P+-k + -q + , P~, -k T ). 

We work in axial gauge n ■ A = for the analysis below, n being a gauge 
vector. In the end-point region with q + — > P + , the s quark has a large 
minus component P~ but a small invariant mass P 2 . Hence, photon vertex 
corrections and self-energy corrections to the s quark in fig. lb give double 
logarithms from the overlap of collinear and soft divergences. In the leading 
regions with loop momentum I parallel to P s and with soft I, the integrand 
associated with the photon vertex correction has the partial expression 

(/P- jk+ /+m b ) lx (f + m B ) 2P x - lx (f-m b )+ / 7a 

[{P-k + iy-m 2 } 2P~l (r + m B ), (1) 

where the factor (/P + m#) is the matrix structure of S, and those terms 



proportional to k have been neglected. The second term in the numerator of 
eq. (0) is suppressed by power (tub — iTib)/ m B- The third term is negligible 
for soft I. For l//P s , it can be easily shown that the trace of 7-matrices 
associated with fig. la, which contains the third term, vanishes. 

Therefore, in the leading regions of / radiated gluons are detached from 
the b quark and collected by an eikonal line described by the Feynman rule 
P\/P-l. The factor 1/P-l is associated with the eikonal propagator, and the 
numerator Pa is assigned to the vertex, where a gluon attaches the eikonal 
line. This observation is consistent with the flavor symmetry in the heavy 
quark effective theory [jDJ. We then absorb this type of radiative correc- 
tions into a jet function J(P S ), which is now factorized out of the process. 
This absorption is reasonable because the soft divergences in fig. lb cancel 
asymptotically as explained later, and the double logarithmic corrections are 
mainly collinear. 

Self-energy corrections to the b quark and loop corrections with real glu- 
ons connecting the two b quarks in fig. lc contain only single soft logarithms, 
and are grouped into the distribution function, which is supposed to be dom- 
inated by soft dynamics. Remaining important corrections from hard gluons 
are absorbed into a hard scattering amplitude H. At last, the factorization 
formula for the spectrum of the decay B — > X s ^ is written as 
dT r 

= J dk + d 2 k T S(k, n)H(P - k, fi)J(P s , fi) , (2) 

which is described graphically in fig. 2. Here E y = g + /v / 2 is the photon 
energy, and fi is the factorization and renormalization scale. 

The Wilson-loop formalism for the decay B — ► X s ^j [[Tj| is summarized 
in fig. 3. Soft gluons radiated from the b quark are collected by the eikonal 
lines along the classical trajectories of the b and s quarks, which form a 
Wilson loop, and are absorbed into S. Since the s quark is almost on-shell in 
the end-point region, the Wilson loop introduces extra collinear singularities 
into the distribution function, which are not consistent with the dominant 
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soft dynamics in the B meson. Note that the light s quark is not replaced 
by an eikonal line in our treatment. Hence, both S and J in [nj contain 
double logarithms, which are treated by RG methods. However, due to 
the presence of double logarithms, relevant anomalous dimensions involved 



in RG equations must be accompanied by large logarithms |L2]], such that 
perturbative calculation is not reliable. 

For the purpose of demonstrating the equivalence of our approach to the 
Wilson-loop one, it is convenient to choose a time-like gauge vector n = P, 
instead of a space- like n as in ||, |, [7|], under which radiated gluons do not 
attach the eikonal lines associated with the b quarks. This choice of n does 
not bring new pinch singularities into S. 

We drop the intrinsic dependence in S, and Fourier transform eq. (0) 
into b space, b being a conjugate variable to fcy, 

dE 1 J "'" (2tt) 5 

It has been shown that integrands for soft corrections are proportional to 
1 — e tlj " b [I], 0], which vanish in the asymptotic region with b — > as stated 
before. 

J depends on P s and n through the ratio (P s ■ n) 2 /n 2 due to the scale 
invariance in n of the gluon propagator in axial gauge, 

To sum up the double logarithms in J, we consider the derivative dJ/dhaP~ 
[Q. 171] , P~ being regarded as a variable from now on. Replacing d/d In P~ by 
d/dn using a chain rule, we have 

d ~ v? d ~ 

- J = -— v «^- J ( 5 ) 



d In P s v ■ n dn a 

with the vector v = (0,1,0). d/dn a operates only on a gluon propagator, 
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and gives |I| 

J_ N ^ = _ J_( N my + N ™^\ , ( 6 ) 
dn a I ■ n 

Combining eqs. and (Bp, we find that a gluon vertex, after contracted 
with I that locates at both ends of the differentiated gluon line, has been 
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replaced by a new vertex, described by the Feynman rule - —, v a . This new 
vertex will be represented by a square. Adding together all diagrams with 



different differentiated gluon lines and using the Ward identity [lq| , this 
square moves to the outermost ends of J. 

The gluon momentum I flowing through the square vertex does not lead 
to collinear divergences because of the denominator n ■ I. The leading regions 
of I are then soft and ultraviolet, in which the subdiagram containing the 
square can be factorized. We then derive a differential equation, 

d 



dlnP: 



- J — 2 [K,{oi s {jj)) + g(a s (/i))] J , (7) 



as shown in fig. 4a, where the eikonal lines disappear due to the choice n = P. 
The factor 2 counts the two ends of J. The functions K, and Q collect the 
soft and ultraviolet divergences in the subdiagram, respectively. Lowest-order 
diagrams of /C are exhibited in fig. 4b, and those of Q are in fig. 4c. 

A straightforward calculation gives K = fig. 4b — 5)C and Q = fig. 4c — 5Q, 
5K, and 5Q = —5JC being the corresponding counterterms. Since /C contains 
only single soft logarithms and Q contains only single ultraviolet logarithms, 
they can be treated by RG methods: 

IL^-K = -X K = -H-^-G ■ (8) 
dn d/j, 

Xfc = /^^7 is the anomalous dimension of /C, whose expression up to two 
loops is given by [p], 

a s „ fa s \ 2 „ /67 tt 2 \ 5 
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with rif the number of quark flavors, and Cp = 4/3 and Ca = 3 the color 
factors. \)c is independent of the gauge vector n, and thus eq. (|J) is identical 
to that derived from a space-like n in J7|. 

There are two P~-dependent scales involved in J: one is P~ itself, which 
serves as the lower limit of /i, and the other comes from the invariant P-P s ~ 
P + Pg ■ Solving eq. (H), we have 



ic + g = ic{a s {p-)) + g 



P+P: 



+ 



— A/c(a a (/x)) 



(10) 



In the second expression we have neglected the initial condition of K, + Q. 
Substituting eq. (|T0|) into (^), we derive 



J{ p ti p s ) b >)") = ex P 



i/6 ^; 

xJ(P/,b,/i) , 



— X K (a s (n)) 



where we have taken the infrared cutoff 1/6 (conjugate to &t) as the lower 
bound of P~ and the largest scale P~ involved in the process as the upper 
bound. The four scales appearing in eq. (pT|) satisfy the ordering P > 



P + Ps > Pg > 1/b > Aq CD , which justifies the neglect of the initial 
condition of K + Q. 

The function J(P+, b,/i) obeys the RG equation, 



Vj(P+,b,v) = -2\jj(P+,b,ti) 



(12) 



with V = /i-^- + f3(g)-§^ and Xj = —a s /ir the quark anomalous dimension 
in axial gauge. The invariant mass of the jet, P s 2 ~ 2P S + P _ « 2P _ /6, if 
assuming that the small P s + is of the same order as 1/b, is a natural initial 
scale of the evolution. The solution to eq. ( |12"D is then given by 

djl 2 



J(P/,b,/i) = exp 



p-/b /i z 



-\j{a s {\x)) 



J(P 



P-/6) . (13) 
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To sum up next-to-leading logarithms, we have to organize the single soft 
logarithms in S. The RG equation for S is written as 



VS(k + , f i) = -X s S(k + ,fi) 



(14) 



where Xs = — (a s /7r)Cp is the anomalous dimension derived from fig. lc with 
the eikonal approximation. We allow S to evolve to the infrared cutoff 1/6, 
and obtain 



S(k + , //) = exp 



—X s (a s (ri)) 

1/6 \X 



S(k + ,l/b). 



(15) 



The initial condition S(k + ,l/b) has a nonperturbative origin arising from 
soft QCD dynamics in the B meson. 

Since the differential decay rate is independent of fi, the RG equation of 
the hard scattering H is written as 



VH(fi) = (2Xj + X s )H(ti) 

which is solved to give 

<p-) 2 dp, 2 



(16) 



H(fM) = exp 



2 Xj(a 8 (fi)) - I ^X s (a s (ii)) 



H(P~ 



We have chosen the largest scale P~ as the upper bound of the integral. 
Combining eqs. (|TTD - (jTT|) , we derive the differential decay rate 



(17) 



dT 1 



d 2 b 



(2tt) 2 

x exp(-s 7 (P~, 1/6)) , 



S(k + , l/b)H(P + - k + } b, P-)J(P+, b, dP-jb) 



(18) 



with the Sudakov exponent 

, / — A/cM/z)) 

1/6 P s JP S - H 

P~ djl^ ,^ , f(P-) 2 dfi 2 



+ / -4x s (a s (fl)) + [ 

Jl/b jJ, JF 



p-/b n' 



\j{a„(ji)) . 



(19) 
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Employing the variable changes P s = yP in the first term, /i = yP in the 
second term and /2 2 = yP~ 2 in the last term, eq. (|T9"| ) reduces to 



dy 

i/bP- y 



2 l^ P ltx K (a s (ji)) + \ s (a s (yP-)) + \j{a s {^yP~)) 

yP- /I 



(20) 



The corresponding expression for the exponent in [15[ is quoted below: 

Ino/N y 



S N 



s/y M dk t 

^ I ~ J- cusp! 

yM k t 



a 8 (kt)) + T(a s {yM)) + 7 (o; s (V£M)) 



(21) 

where sn describes the evolution of the iV-th moment of the differential decay 
rate, with M the b quark mass, the constant n Q = e _7E , 7^ being the Euler 

-(a s /ii)CF and 7 = —a s /7[ the anomalous 



constant, and T 



cusp 



Ak, r 



dimensions. The two-loop expressions for r cusp and T can be found in [12, IS 



Comparing eq. (p0|) to (|21|), we easily identify the equality of the anoma- 
lous dimensions: 



A 



K, 



CUSp 7 



A.. 



7 



(22) 



The two expressions s 7 and s n are basically the same except the lower bounds 
of y. The scale 1/b comes from the inclusion of transverse degrees of freedom 
in our formalism, and the iV dependence is due to the moment analysis in 



3. The Semileptonic Decay B — ► X u tv 

The analysis of the inclusive semileptonic decay B —>■ X u lv is similar to 
that of the radiative decay B — > Xg'j. Large double logarithmic corrections 
are also observed near the high end of the charged lepton spectrum. We 
choose the B meson momentum P and the light parton momentum k as in 
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Section 2. The lepton momentum Pi has only a large plus component P^~. 
The neutrino carries the momentum P u , Pj = 0. The u quark, emitted by 
the b quark through the transition u / y fl b, has the momentum P u = P — k — q 
with the lepton pair momentum q = Pe + P v = (P/ + P+, P~, P v t)- When 
P/ reaches its maximum P + , implying P+ m and ~P v t ~ 0, the u quark 
has a small invariant mass, and moves fast in the minus direction with P~ = 
P~ — P~ . The situation is then similar to the end-point region of the radiative 
decay B — > X s j. 

Repeating the resummation procedures in Section 2, we derive the cor- 
responding Sudakov exponent s^, which is similar to s 7 but with the largest 
scale P~ replaced by (1 — x v )P~ , x u = P~ / P~ . se is given by 



Ji/b P u Jp- fj. 

+ / , -^XsMfi)) + -^-\j(a 8 (jA)) 

Jl/b jJL J(l-x v )P-/b jJL z 



Jl/bP- y JyP- U 



"Vvp- dp 

lyP- p 



+ / 2L\j{a.(y/yP-)) . (23) 

J(\-x v )jbP- y 

Note the change in the ^-integration limits compared to eq. (|20|) . 



4. Conclusion 

In this letter we have shown in detail that the resummation technique 
is indeed equivalent to the conventional Wilson-loop formalism. The two 
approaches lead to the same resummation results in the end-point region 
of inclusive B meson decays, except that the large logarithm appears as 
ln&P - in our expression and as In TV in ||15|| . This equivalence provides a 
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justification of our PQCD analysis of B meson decays including Sudakov 
effects in previous works [|, |, 0]. 

We have neglected the initial condition of K, + Q in demonstrating the 
equivalence, which in fact contributes to next-to-leading logarithms. How- 
ever, we believe that physical predictions should be insensitive to the dis- 
tinction in next-to-leading logarithms. It is known that a substitution of the 
scale P~ by another one CP~, with C a constant of order 1, is allowed in 
resummation |L], [18] . This constant C in the integral involving Xfc then makes 
a difference at the level of next-to-leading logarithms. Hence, if the initial 
condition of K + Q is included, predictions will be modified only slightly. 
Certainly, this issue needs more careful study. Other interesting subjects, 
such as the gauge dependence of resummation and the application of our 
technique to deep inelastic scatterings and Drell-Yan productions, will be 



discussed in a seperate work pO 
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Figure Captions 



Fig. 1. (a) Basic factorization for the decay B — > X s 7, and (b) and (c) 
0(a s ) radiative corrections. 

Fig. 2. Factorization for the decay B — ► Xgj in our formalism for a general 
gauge vector n. 

Fig. 3. Factorization for the decay B — > X s ^ in the Wilson-loop formalism. 

Fig. 4. (a) Graphic representation of eq. (|7|) for a gauge vector n = P, and 
0(a s ) diagrams for (b) the function /CJ and for (c) the function QJ. 
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